The Realization of Artificial Kondo Lattices in Nanostructured Arrays 
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The interplay of magnetic energies in a Kondo lattice is the underlying physics of a heavy fermion 
system. Creating an artificial Kondo lattice system by localizing the moments in an ordered metallic 
array provides a prototype system to tune and study the energetic interplay while avoiding the com- 
plications introduced by random alloying of the material. In this article, we create a Kondo lattice 
system by fabricating a hexagonally ordered nanostructured array using niobium as the host metal 
and cobalt as the magnetic constituent. Electrical transport measurements and magnetoresistivity 
measurements of these artificial lattices show that the competing exchange coupling properties can 
be easily tuned by controlling the impurity percentage. These artificial Kondo lattice systems enable 
the exploration of an artificial superconductor which should lead to a deep understanding of the 
role of magnetism in unconventional superconductors. 



The role of magnetism in unconventional 
superconductors-cupratesj 1 -' heavy fermions^ and 
the recently discovered Fe-based pnictidePI -is a current 
area of intense debate for physicists and material 
scientists alike. Typically, the energy scale of magnetic 
excitations and the superconducting order parameter 
have a one-to-one correspondence. 4 5 In heavy fermion 
compounds, this magnetic energetics is reflected via 
the interplay of the localized Kondo interactional 
the screening interaction between a single magnetic 
atom and surrounding conduction electronsand the 
long-range magnetic order. 8 Heavy fermion compounds 
of intermetallic rare earth materials have the natural 
tendency to create the Kondo interaction like phenomena 
between a chemically ordered lattice of local moments 
(/-band) and the compounds conduction electrons (s 
and p-bands)P 

The formation of Kondo lattice in rare earth inter- 
metallic compounds (of cerium, ytterbium, uranium) 
lead to many surprising phenomena which includes the 
observation of unconventional superconductivity (for ex- 
ample UPd2Al3, CeRhlns/^^an d qua ntum critical phe- 
nomena (for example YbRh2Si2)! 13 * 14 * Quantum critical 
phenomena (QCP) are accompanied by a second order 
quantum phase transition between a long range mag- 
netic order regime, favored by Ruderman-Kittel-Kasuya- 
Yoshida (RKKY) interaction which results from an indi- 
rect exchange interaction between local moments medi- 
ated by conductio n elec trons, and the Kondo screening 
regime at T = KpEHThe competing energy scales in- 
volved in QCP and their scaling with temperature is ar- 
gued to be directly related to the superconducting order 
parameter, gap energy A, in some of these compounds 
and therefore the magnetic energy interplay is considered 
at the core of the occurrence of superconductivity.^ The 
crystalline nature of a Kondo lattice provides an ideal 
setup for study. In this article we describe the creation 
of artificial magnetic crystals in nanostructured arrays. 
The hexagonal ordering of magnetic impurities (cobalt) 
in the host metal (thin niobium film) automatically cre- 
ates a magnetic crystal. Controlling the percentage of 




FIG. 1: (color online) Schematic description of artificial mag- 
netic crystals, (a) As we see in this figure, localized magnetic 
ions (dark grey) embedded in a conducting film form a hexag- 
onal lattice pattern and thus create an artificial magnetic crys- 
tal. If the amount of magnetic ions at each scattering site is 
small then each magnetic scattering site is totally screened 
by conduction electrons and therefore forms a dilute Kondo 
lattice system, (b) When the percentage of magnetic ions 
at each scattering site is increased with respect to the con- 
ducting film then indirect exchange interaction, mediated by 
conduction electrons, between localized magnetic ions plays 
a significant role. In figures (a) and (b), the percentages of 
magnetic ions at scattering sites are reflected by the (Red) 
arrows weight. 

magnetic impurities embedded in each site in the host 
metal acts as the tuning parameter to allow for the de- 
tailed investigation of the energetic interplay between 
Kondo and RKKY-type interactions (Fig. 1). Exchange 
coupling constants were estimated using electrical and 
magnetoresistance measurements of the artificially fab- 
ricated Kondo lattices. As shown in the schematics in 
Fig. 1, our technique leads to the implantation of Co 
particles inside the Nb film which is a different physical 
system from magnetic dot arrays on top of a substrate. 
In the latter case, conduction electrons do not interact 
with the magnetic impuritys moment directly. 



The first step in the fabrication process of artificial 
Kondo lattices involves the development of hexagonal 
nanoporous copolymer templates from a self-assembled 
diblock copolymer film of thickness 36 nm and aver- 
age pore diameter of 13 nm on top of a 10 nm thick 
conducting substrate (Nb). Lattice separation between 
pores was 28 nm. Details about the method of devel- 
oping the nanoporous copolymer template can be found 
elsewhereP^l The angular deposition of the desired mate- 
rial onto the walls of the nanopores was followed after this 
procedure. During the angular deposition process, the 
substrate (Nb film with nanoporous polymer template 
on top of a silicon wafer) was uniformly rotated about 
the axis perpendicular to its plane. Angular deposition 
depends on a critical deposition angle C = tan _1 (D/h), 
where D and h are the diameter and height of a nanopore, 
respectively. In our experiment, a deposition angle of 6 = 
25 degree (> 6 C = 20 degree) was chosen to avoid the de- 
position of Co at the bottom of the pores. The thickness 
of deposited Co material on the walls of nanopores was 
in the range of 4-5 nm (based on the calibrated quartz- 
crystal microbalance reading). This technique of angular 
deposition of Co onto the wall of the nanopores results 
in a multiply-connected geometry as the top of the film 
(polymer template) is always exposed to the evaporant, 
as shown in Fig. 2a and 2b. After material (cobalt) 
deposition, the calibrated parallel beam ion-beam etch- 
ing technique was used to forcefully knock the desired 
amount of Co particles out of the nanopores walls and 
embed into the Nb substrate in an ultra-high vacuum 
chamber filled with argon gas at 10 -5 torr. For this pur- 
pose, the sample was tilted by 20° during the ion-beam 
etching and the ion beam was accelerated at 2000 eV. 
A schematic description of this process is shown in Fig. 
2c-d. By controlling the etching rate and time, we were 
able to control the amount of Co particles embedded in 
a hexagonal order in Nb film. After etching, the sample 
was rinsed with toluene for 6 hrs to remove the remain- 
ing polymer template encapsulated with thin cobalt film, 
leaving only the Nb film with magnetic impurities. Fig. 
2e-f show TEM (transmission electron microscope) im- 
ages of artificial lattices resulting after 6 and 10 minutes 
of ion-beam etching. We clearly see in these figures that 
a certain percentage of Co particles are embedded in the 
host substrate and form a hexagonal matrix. For conve- 
nience, samples resulting after different etching times are 
named Nb-diblock (etching time = min, equivalent to 
Fig. 2a), NCI (etching time = 2 min), NC2 (etching time 
= 4 min), NC3 (etching time = 6 min) and NC4 (etch- 
ing time = 10 min). Based on the calibrated ion-beam 
etching rate, it is roughly estimated that samples NCI, 
NC2, NC3 and NC4 have no more than 0.03%, 0.06%, 
0.1% and 0.2 % of Co particles with respect to Nb, re- 
spectively. In terms of the number of Co particles per 
site, most dilute sample NCI has about 25 Co particles 
per site. The uncertainty in the determination of the ex- 
act percentage of magnetic impurities embedded in the 
conducting substrate does not affect our conclusion. 
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FIG. 2: (color online) Fabrication scheme and electron micro- 
scope images of artificial Kondo lattices, (a) & (b) Schematic 
views from top and side of multiply-connected hexagonally or- 
dered diblock copolymer porous template encapsulated with 
cobalt film are shown in these figures. In this figure, or- 
ange color posts are diblock copolymer walls, white spaces are 
empty pores and blue capping is cobalt film, resulting from 
the angular deposition of cobalt on the wall of the pores, (c) 
Calibrated ion-beam etching technique is used to knock out 
desired amounts of cobalt particles from the encapsulating 
cobalt layer and forcefully implant them in the conducting 
substrate, (d) This schematic figure shows the top view of 
conducting film with cobalt particles embedded in hexagonal 
matrix after the ion-beam etching, (e) TEM image of sample 
NB3 (6 min. etching) is shown in this figure. We clearly see 
local hexagonal ordering of magnetic impurities in this figure, 
(f) TEM image of a large array of 10 min ion-beam etched 
film (NB4). 

Electrical transport measurement data in the temper- 
ature range 4-300 K are shown below in Fig. 3a. Below 4 
K, the substrate (Nb film with diblock polymer template 
only) becomes superconducting and no useful informa- 
tion regarding the magnetic properties of the system can 
be extracted.^ Therefore we will focus on electrical mea- 
surements above 4 K only. Electrical resistivity data in 
Fig. 3a were normalized with respect to the value at T = 
290 K to get a better comparison between different sam- 
ples. Comparing the electrical transport data of all five 
samples, we observe different variations of the resistivity 
with temperature. The resistivity of the undoped Nb- 
diblock sample has been reported before and we will not 
discuss it here. 18 Sample NCI exhibits a deep Kondo-like 
resistance minimum around 70 K and saturates below 30 
K as the measurement temperature is further reduced. 
The saturation of resistivity below T = 30 K suggests 
the onset of Kondo screening and formation of the Kondo 
lattice in this system.^ This will be discussed in greater 
detail in the following paragraphs. In the case of NC2, 
the resistivity first increases, peaks around 200 K and 
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then decreases before increasing again below ~ 40 K. The 
presence of the resistance minimum at T~30 K suggests 
that the Kondo effect is still present, perhaps weakened, 
in this system. The broad maximum in resistance around 
200 K is difficult to analyze due to possible contributions 
from many different effects such as thermal fluctuations 
and spin-diffusion. Nonetheless it is reasonable to argue 
that the resistance maximum at c± 200 K corresponds to 
local spin correlations which exceed thermal disorder and 
begin to form magnetic impurity clusters via an indirect 
exchange interaction!^ No Kondo like resistance mini- 
mum is observed in sample NC3; instead the resistivity 
starts increasing as the measurement temperature is re- 
duced untill T~100 K. Below 100 K, the resistivity drops 
significantly before saturating below T~30 K. Similar be- 
havior of resistivity variation has been observed before in 
Kondo lattice systems and the onset temperature is often 
associated with the characteristic temperature indicating 
the transition to long range magnetic order mediated by 
the RKKY-type interaction! 21 * 22 * Further increase in the 
percentage of magnetic impurities as a result of longer 
ion-beam etching leads to larger local moments. In that 
case, direct magnetic interaction between localized mag- 
netic impurities starts playing a role. This sets up a new 
regime marked by the absence of Kondo like interaction 
and the presence of direct magnetic interaction between 
localized impurities, as observed in sample NC4 where 
an increase in resistivity at high temperature is followed 
by a gradual decrease as the temperature is lowered. To 
further elucidate the properties of magnetic crystals and 
learn about the interplay of Kondo interaction vs. RKKY 
interaction, the following paragraphs in this article will 
primarily focus on the resistive and magnetoresistive be- 
havior of samples NCI, NC2 and NC3 only. 

The Kondo interaction leads to the formation of com- 
posite fermions with spin quantum number S=1/2P^ 
Theoretically it has been shown by Costi et al. 
that Wilsons numerical renormalization group (NRG) 
method 2 ^ can be effectively used to extract the character- 
istic Kondo temperature from transport measur emen ts of 
dilute magnetically impure conducting systems! 24 * 25 * The 
application of the NRG method to the Kondo model char- 
acterizes the physical properties of a Kondo system in 
which the impurity moment weakly couples to the con- 
duction electrons. In Fig. 3b, we have fit the low temper- 
ature resistivity data of sample NCI using the NRG solu- 
tion for S = 1/2. As we see in this figure, the NRG solu- 
tion for S = 1 /2 gives a reasonable fit to the experimental 
data with = 56 K. To further verify the Kondo transi- 
tion temperature, we have also fit the experimental data 
using the Kondo expression: p~-ln(T/Tx) in the tem- 
perature range of T = 70 K to T = 30 K. The obtained fit 
of Tr-j 56.7 K, is consistent with the NRG fitting value. 
These fits suggest that sample NCI forms a Kondo lat- 
tice system with Kondo screening temperature of ^ 56 K. 
In the present analysis of experimental data, we ignore 
the weak-localization (WL) and electron-electron inter- 
action (EEI) effects. In general, WL and EEI effects are 
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FIG. 3: (color online) Resistivity measurement analysis of 
magnetic crystals, (a) Normalized resistivity data of different 
magnetic crystals are shown in this figure along with the resis- 
tivity data of niobium film with diblock copolymer template 
only (reference sample). As the magnetic impurity percent- 
age is increased at each scattering site, the physical behavior 
changes from Kondo screening (NCI) to RKKY-type interac- 
tion (NC3). (b) Resistivity data of sample NCI is fitted with 
NRG soln. for spin-1/2 and Kondo expression. Obtained 
fitting parameters are reasonably consistent with each other 
(see text). Also shown in this figure is the log- normal behavior 
(Red curve) exhibited by NCI over broad temperature range, 
(c) Low temperature data of sample NC2, fitted with Kondo 
expression, is shown in this figure, (d) Resistivity data and 
resistivity derivative data of sample NC3 are shown in this 
figure. A clear peak is observed in resistivity derivative data 
at T = 51 K while a large drop in resistivity starts occurring 
around T~70 K. 

more dominant in quasi one-dimensional systems where 
the dimensional length scale is smaller than the mean 
free path for conduction electrons, which is not the case 
in this work. Interestingly, the NCI resistivity curve is 
also well fitted with log-normal behavior over a wide tem- 
perature range, as shown by the red curve in Fig. 3b. 
The obtained fitting value of = 59.2 K in this case 
is slightly larger than other fitting results for this sys- 
tem. This unusual log-normal behavior of resistivity as a 
function of temperature is observed for the first time in 
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Kondo systems and suggests that Kondo screening and 
logarithmic normalization of resistivity are inter-related. 

We have also estimated the exchange interaction be- 
tween the Co magnetic impurities and conduction elec- 
trons of Nb using the Kondo expression KbTk — exp(- 
1/ Jn)^ where J is the absolute value of exchange interac- 
tion constant and n is the density of states of conduction 
electrons. For n=5.5xl0 28 m _3 J^the estimated value of 
J is found to be -9.7 meV. The estimated exchange con- 
stant for sample NCI is of the same order as observed 
in other dilute Kondo lattice systems!^ We have also 
fit the low temperature experimental data of NC2 using 
the logarithmic Kondo equation, as shown in Fig. 3c. 
The best fit is obtained for Tk — 24 K which corre- 
sponds to J = - 6.3 meV. It is noteworthy to mention 
that each hexagonal magnetic impurity location in NC2 
has a higher percentage of Co as compared to NCI. That, 
possibly, prohibits the application of Kondo behavior to 
this system as the Kondo formalism is more applicable 
to dilute impurity cases. Increasing the concentration of 
magnetic impurities also leads to the onset of indirect 
exchange interaction between impurities. 

Sample NC3 does not exhibit Kondo-type minima in 
resistivity measurements; rather a large drop in resis- 
tivity is observed around 75 K. We have also plotted 
dp/dT vs. T in Fig. 3d. A clear peak is observed in 
differential resistivity data at T = 51 K which quantifies 
the point of inflection of the resistivity curve. Similar 
behavior in resistivity measurements has been observed 
before in concentrated Kondo lattice systems and were 
associated with the RKKY-type indirect exchange in- 
teraction between magnetic impurities mediated by con- 
duction electrons. 28 Since RKKY-type interactions be- 
tween magnetic ions leads to long range magnetic order, 
many researchers argue that the resistivity maxima rep- 
resent the magnetic transition temperature, T^. 21 Since 
we are interested in estimating the exchange interaction 
strength of RKKY-type interaction for sample NC3, we 
chose the point of inflection of the resistivity curve, Fig. 
3d, as the characteristic temperature, Trkky = 51 K, 
below which conduction electron mediated RKKY-type 
interaction plays a significant role. RKKY exchange cou- 
pling strength, Jrkky, is estimated using the expression 
KbTrkky — J 2 nW^ The estimated magnitude of the 
exchange coupling constant Jrkky is approximately 2.1 
meV. Thus we see that the exchange coupling strength 
changes as the amount of magnetic impurities is varied 
from dilute (2 min. etching) to concentrated (6 min etch- 
ing) in the hexagonally localized matrix. 

To further understand the magnetic properties of 
Kondo lattices, we have performed magnetoresistance 
measurements for both in-plane and perpendicular ap- 
plication of magnetic fields with respect to the plane of 
the substrate. In Fig. 4a-c, the magnetoresistance data 
of NCI, NC2 and NC3 are plotted as a function of tem- 
perature. Magnetoresistance (MR) is defined as MR = 
[p(H) - p(0)]/p(0). As we can see in Fig. 4a, the applica- 
tion of an external magnetic field up to 2 Tesla on sample 
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FIG. 4: (color online) Magnetoresistance (MR) measurements 
of Kondo lattices, (a) Magnetoresistance data of sample NCI 
for both in-plane and perpendicular applications of external 
magnetic field are plotted in this figure. No clear difference 
is observed between measured data at different applied fields, 
indicating complete Kondo screening of localized moments, 
(b) In the case of sample NC2, MR data for perpendicular 
application of magnetic fields lie higher compared to in-plane 
data of same field values. It suggests that Kondo screening 
is weakened and the direction perpendicular to the substrate 
develops as the easy axis for this system, (c) In this figure, 
we have plotted the MR data of sample NC3. Significantly 
higher MR percentages are observed at low temperature in 
this case. We also observe a deep minimum in MR indicating 
the onset of coherence and heavy fermion state (see text), (d) 
We have plotted the change in resistance, AR, as a function of 
perpendicular magnetic field at T = 7 K in this figure. Solid 
line is a guide to the eye. Observed peaks in the AR data 
are attributed to the cyclotron radius of A-B rings circling 
different number of scattering centers (as shown in the inset). 

NCI does not change the MR substantially. Also, no ob- 
servable differences between MR data for in-plane and 
perpendicular application of magnetic fields are found. 
This simply indicates that the magnetic impurity spins 
are completely screened by surrounding conduction elec- 
trons and hence remain unaffected by the application of 
an external magnetic field. On the other hand, the mag- 
netoresistance data of sample NC2 (Fig. 4b) exhibit dif- 
ferent responses for magnetic field applied in-plane and 
perpendicular to the substrate. MR value at a partic- 
ular temperature is larger at higher magnetic field and 
overall, magnetoresistance curves for perpendicular field 
application lie higher compared to MR curves for in-plane 
field application. The perpendicular direction to the sub- 
strate, in this case, evolves as the easy axis. 

In the MR measurements of NC3, two interesting 
behaviors are observed: the formation of MR minima 
around ~ 80 K which suggests the onset of coherence 
and heavy fermion state^ and the fluctuation of mag- 
netoresistance as a function of perpendicular magnetic 
field. MR value for perpendicular application of field in 
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this case is as high as ~20 % at low temperature. In Fig. 
4d, we have plotted the change in resistance AR = R(H)- 
R(0) as a function of perpendicular applied field at T — 7 
K. We clearly see the oscillatory nature of AR as a func- 
tion of field in this figure. In general, localized magnetic 
impurities act as the scattering centers for conduction 
electrons which lead to interference effects between con- 
duction electrons when the phase coherence of electrons 
exceeds the inter-impurity distance.^ If the magnetic im- 
purities are arranged in a regular fashion then conduction 
electrons form Aharonov-Bohm (A-B) rings encircling a 
discrete number of impurities. 32 Radius (tab) of A-B 
rings are given by AB± ~ h/ '(irer ab 2 )- In Figure 4d, we 
clearly see a strong peak in magnetoresistance data at 
0.15 T and a relatively weaker peak at H~l T. Using 
the above expression for A-B rings radius, these magnetic 
field values correspond to A-B rings of 300 nm and 120 
nm radiuses respectively. These cyclotron radiuses cor- 
respond to electrons encircling 64 and 24 magnetic scat- 
tering centers respectively, as shown in the schematic de- 
scription in the inset of Figure 4d. Recently Siegert et al. 
identified similar experimental observation in a semicon- 
ductor heterostructure with RKKY-type interaction by 
correlating the oscillatory behavior of indirect exchange 
constant, J, with the A-B radius tab via range function 
over a large distance (^500 nm)P 

The energetic interplay between the small exchange 
energy Kondo screening regime and the large exchange 
energy RKKY-type regime can also be easily tuned by 
changing the lattice separation between magnetic impu- 
rities. The nature of RKKY indirect exchange interac- 
tion varies between ferromagnetic and antiferromagnetic 
interactions with distance! 3 -^ Thus the combination of the 
change in magnetic impurities percentage with the vari- 
ation in lattice separation can be used to tune both the 
strength and nature of energetic interplay. Experiments 
of artificial Kondo lattices were also repeated with alu- 
minum as the host metal and similar resistivity patterns 
were observed. However, the estimated exchange con- 
stants in different regimes (Kondo screening and RKKY- 
type) were found to be relatively larger. As the thick- 
ness of the metallic host was increased to 25 nm, a sharp 
downturn in the electrical resistance data was observed 
below the resistance maxima. This suggests that the en- 



ergetic interplay can be further explored by choosing dif- 
ferent combinations of host metal and magnetic impuri- 
ties. It will be interesting to explore the scaling behavior 
in these artificial Kondo lattice systems, in comparison 
to those recently reported by Yi-Feng et al. in inter- 
metallic rare earth compounds.^ The scaling of Kondo 
lattices involve a characteristic temperature, T*, defined 
by the intersite coupling energy! 3 -^ Precise estimation of 
T* in these magnetic crystals would require magnetic sus- 
ceptibility and heat capacity measurements. 

We now discuss the important implication of artificial 
Kondo lattice creation. Many Kondo lattice compounds 
exhibit the competing dual phenomena of unconven tional 
superconductivity and magnetic energy interplay! 9 * 10 * 12 * 
However, the observation of superconductivity is not uni- 
versal in all Kondo lattice systems. This prohibits any 
strong theoretical formulation combining magnetic en- 
ergy interplay and superconductivity. Here we propose a 
reverse approach using artificial Kondo lattice systems. 
We have shown that the energetic interplay of the com- 
peting magnetic interactions occurs in the range of -9.7 
meV to 2.1 meV in the presented artificial Kondo lat- 
tices, similar to that in the bulk Kondo lattice systems 
exhibiting unconventional superconductivity. If the mag- 
netic energy interplay and superconductivity are directly 
related and the magnetic excitations provide a mediat- 
ing mechanism for Cooper pair formation, essential for 
superconductivity, then it would be possible to create an 
artificial crystal with superconducting properties. Ad- 
ditional study of the artificial magnetic crystals using 
different metallic host with higher electronic density of 
states and different magnetic constituents are necessary 
for this purpose. Discovery of superconductivity in an ar- 
tificial Kondo lattice system will not only firmly establish 
the role of magnetism in unconventional superconductors 
but will also provide an important tool to manipulate the 
physical properties of the materials, which will be crucial 
for the future development of higher temperature super- 
conductors. 
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